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I . INTRODUCTION 

In  this  report  we  describe  progress  made  in  three  different  areas 
of  research.  In  all  cases,  our  two  previous  technical  reports  [1,2] 
are  recommended  background,  and  we  shall  freely  call  upon  results 
derived  in  them.  First,  a new  image  quality  factor  Q is  defined,  and 
a variety  of  curves  of  Q vs.  the  average  number  of  image  photoevents, 

N , are  given  for  both  precompensated  and  postcompensated  imagery. 

Second,  a comparison  of  the  photon  noi se- 1 im i ted  performance  of  amplitude 
i nterf eromet ry  and  speckle  interferometry  is  presented.  Lastly,  we 
describe  the  current  status  of  our  formulation  of  the  least-squares 
filtering  problem  when  the  filter  is  allowed  to  be  space-variant  or 
non  1 i near . 


I I . A NEW  QUALITY  FACTOR  Q 

In  previous  reports  we  have  used  a quality  factor  Q defined  by 


Q 


E 


-1 


(0 


where 

N = the  mean  number  of  photoevents  detected  in  the  entire  image; 

*t* 

S = the  transfer  function  of  the  ideal  system,  which  is  taken  to 
be  that  of  a d i f f ract ion- 1 imi ted  telescope  with  a circular 
aperture; 

and  S2y  = spatial  frequencies  measured  in  cycles  per  radian 
of  arc  in  the  sky; 

'here  and  throughout  this  report,  quantities  with  a " over  them  are 
normal i-ed  to  have  value  unity  at  the  origin.  j > 





r-  'T , 1^ 


/N 

B «=  total  blur  OTF 

iJIq  = the  normalized  average  energy  spectral  density  associated 
with  the  object  ensemble;  and 

E = the  mean-squared  error  associated  with  the  restored  image. 

The  mean-squared  error  was  stated  explicitly  as  being  [2] 

f?  % 

£ = . . dS^yd^Y  (2) 

U ' »0 

Note  that  the  numerator  of  the  fraction  in  expression  (1)  is 
simply  the  ideal  signal  power.  It  is  necessary  to  subtract  unity  from 
the  ratio  due  to  the  fact  that  in  the  limit  where  N 0 , the  denominator 
E approaches  the  numerator,  whereas  we  desire  a quality  measure  that 
goes  to  zero  as  N goes  to  zero. 

Previous  expressions  for  the  transfer  function  of  the  restoration 
filter  have  been  presented  in  normalized  form  [1,  p.l6] 


H t ^ 


(l+N)  Sb4^ 
1 + n|b|^$q 


(3) 


In  order  to  arrive  at  a more  satisfying  definition  of  the  quality  factor 
Q v.'e  have  found  it  convenient  to  deal  with  un-norma  1 i zed  quantities  in 
some  cases.  The  form  of  the  un-norma 1 i zed  optimum  transfer  function  can 
be  shown  to  be 


kNSB*4) 

H(f2y,Sly)  = r- o— 

^ ^ I + n1b|^<i>q 

with  k a constant  given  explicitly  by 


(M 


(5) 


2 - 


where 

h = Planck's  constant 
V = the  mean  optical  frequency 
n = the  quantum  efficiency 
T = the  image  integration  time 

* the  volume  under  the  blur  spread  function  b(x,y),  i.e., 


b 


0 


oo 

b (x,y)dxdy 

— oo 


(6) 


When  the  constant  k is  included  in  the  expression  for  the  optimum 
transfer  function  H , the  expression  for  the  mean-squared  error  becomes 


c 


■7  ip 

■1  ' * *0 


(7) 


This  expression  for  c serves  as  a useful  and  meaningful  representation 
of  the  "noise"  associated  with  the  restored  image,  where  the  term  "noise" 
is  used  here  in  a general  sense  to  include  both  random  fluctuations 
introduced  by  the  random  positions  and  random  numbers  of  the  photoevents, 
and  residual  blur  in  the  restored  image. 

The  numerator  or  "signal"  component  of  the  new  expression  for  Q 
is  also  modified  to  be  more  reasonable.  Rather  than  using  the  "ideal" 
signal  power  in  the  image,  we  use  instead  the  actual  signal  power  in 
the  restored  image.  Thus  the  numerator  of  the  expression  for  Q becomes 


S 


“OO 

where  H is  the  un-norma I i zed  restoration  transfer  function, 
quality  factor  Q thus  becomes 


(8) 
The  new 


- 3 - 


i-iWuniMiHiy jii 


nmiijft;^ipiwwi^igi{W^^inwi^i^wiuiiijt]iiwiiiiPPW"in  lumi.iiuiiunii.i  i.immiiuM. 


(N)^ 


BHT  -1>Q  dU^dil^ 
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(n)2 


is|2 ; 


(9) 


1+N  B $, 


Substituting  Eq.(^)  for  H , we  obtain 


Q 


(N)^ 


1B|^|S|^  1] 

/N  /»• 

(i+N|Br  ^q) 


2 


2 ^ 
|S|  ^ 


1+N  B 

-00  f 


dS^^  d^Y 


(10) 


The  new  quantity  parameter  defined  above  is  more  reasonable  in  its 
behavior  than  the  parameter  used  previously.  Its  most  desirable 
properties  are  that  Q ^ 0 as  N +0  and  Q ultimately  increases  in 
proportion  to  N as  N grows  arbitrarily  large.  This  latter  property 
is  expected  for  any  measure  of  the  mean-squared  s i gna 1 -to-no i se  ratio 
of  a Poisson-related  variable. 

We  now  present  some  specific  curves  of  the  new  Q vs.  for 
several  cases  of  specific  interest.  The  cases  treated  here  are  precisely 
the  same  ones  treated  in  our  last  technical  report  [2],  and  the  reader 
may  wish  to  consult  that  reference  for  a more  detailed  discussion  of  the 
assumptions  behind  each  set  of  curves.  In  all  cases  it  is  assumed  that 
the  object  of  interest  is  an  ideal  point  source  (4)^  = 1),  Figure  1 
shows  plots  of  lOgjQQ  vs.  logj^N  for  three  different  cases  of  post- 
compensation.  In  part  (a)  rt  is  assumed  that  a long-t ime-averaged  image 
is  recorded.  Various  values  of  the  atmospheric  coherence  diameter  r^ 
are  shown,  the  case  = “ corresponding  to  a perfect  diffraction- 


- - 
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limited  image,  with  only  photon  noise  limiting  the  restoration.  In  part 
(b) , it  is  assumed  that  a perfect  tilt-removing  mirror  is  used  as  a pre- 
compensator, and  that  "far  field"  atmospheric  propagation  conditions  hold. 
Finally,  in  part  (c) , tilt  removal  is  again  assumed,  but  "near  field" 
propagation  conditions  hold.  Note  the  break  point  in  the  r^  = 20  cm 
curve,  after  which  Q increases  in  proportion  to  N.  All  curves  have 
such  a break  point  for  large  enough  N , but  only  for  the  r^  = 20  cm 
curve  does  this  point  fall  within  the  range  of  N shown  here. 

For  comparison  purposes,  we  show  in  Figure  2 the  corresponding 
curves  of  Q vs.  N when  ri£  post-compensation  is  applied.  Comparison  of 
Figs.  1 and  2 thus  provides  an  indication  of  how  much  gain  in  quality 
factor  is  provided  by  postcompensf t ion  at  each  level  of  N. 

In  presenting  the  results  for  precompensation  alone  and  for  combined 

pre-  and  postcompensation,  the  theoretical  development  oresented  in 

reference  [2]  has  been  followed.  Specifically,  a shearing  interferometer 

with  317  subapertures  has  been  assumed.  The  ratio  of  image  integration 

time  T to  wavefront  sensor  integration  time  t has  been  taken  to  be 
/* 

10  , and  errors  introduced  by  the  temporal  dynamics  of  the  atmosphere 

have  not  been  included.  The  parameter  3 represents  the  fraction  of  photons 
sent  to  the  wavefront  sensor,  while  I-3  represents  the  fraction  sent  to 
the  image  detection  plane. 

The  curves  of  Fig.  3 show  'ogjQQ  vs.  'o9|qN  for  several  cases 
of  precompensation  alone  (i.e.,  no  postcompensation  filtering  used). 

Part  (a)  assumes  that  3 = O.9  and  shows  the  dependence  of  Q on  N 
for  several  values  of  the  atmospheric  coherence  diameter  r^.  Part  (b) 
assumes  that  r^  * 10  cm  and  shows  the  corresponding  dependence  for  three 
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to 
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Fig.  3:  Q vs.  N with  precompensation  alone: 

(a)  8*0.9,  several  values  of  r^; 

(b)  rQ*10cm,  several  values  of  8; 

(c)  rQ=“>  , several  values  of  8- 


S 00 


values  of  6 (0.3i  0.9.  0.99).  Part  (c)  assumes  that 

(no  atmosphere)  and  again  shows  the  dependence  of  Q on  N for  the 
same  three  values  of  3.  Note  in  the  curves  of  part  (b) , in  the  absence 
of  any  post-compensation  the  curves  appear  to  saturate  at  a finite 
value  of  Q , a property  attributable  to  the  wavefront  fitting  errors, 
which  do  not  depend  on  N.  Such  errors  are  not  present  when  > 

hence  the  curves  of  part  (c)  continue  to  increase  with  N.  Note  also 
that  for  the  curves  of  part  (b),  the  large  values  of  3 (0.9  and 

0.99)  achieve  the  best  values  of  Q for  a given  N , but  when  no 

atmosphere  is  present  (part  (c)),  better  image  quality  is  ultimately 
achieved  (for  large  N)  if  less  light  is  diverted  to  the  wavefront 
sensor.  Since  in  the  latter  case  the  wavefront  sensor  serves  no  useful 
purpose,  these  results  make  good  sense. 

Figure  4 shows  several  results  of  interest  for  combined  precompensation 
and  postc  , sat  ion.  In  part  (a)  of  the  figure,  r^  is  taken  to  be  10  cm, 

and  the  dep  dence  of  Q on  N is  shown  for  several  values  of  3.  The 
curves  do  not  saturate  in  this  case  at  large  N,  due  to  the  fact  that  post- 
compensation removes  the  effects  of  wavefront  fitting  errors.  In  part 
(b)  of  the  figure,  log^gQ  is  plotted  vs.  the  splitting  ratio  3 for 
several  fixed  values  of  N.  At  high  light  levels  (large  N) , an  extremely 

broad  maximum  is  found,  so  the  exact  value  of  3 is  not  of  great  conse- 

quence. At  low  light  levels  (smaller  N) , a narrower  maximum  develops 
in  the  region  of  large  3 (i.e.,  most  of  the  light  sent  to  the  wavefront 

sensor).  Finaily,  in  part  (c),  we  have  log|Q<l  vs.  for  a fixed 

splitting  ratio  (3  = 0.9)  and  for  several  values  of  r^.  The  curves 
for  r^  = “>  and  r^  = 20  cm  are  almost  indistinguishable,  while  Q is 
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noticeably  reduced  by  smaller  r.  , 

U 

fitting  errors. 


a consequence  of  the  increased 


Lastly,  in  Figure  5 we  plot  several  curvp-.  of  logj^Q  vs.  log^^N 
for  comparison  purposes.  In  all  cases  r^  o taken  to  be  10  cm.  Curve 
(a)  corresponds  to  the  case  of  both  precompensation  and  postcompensation, 
with  e “ 0.9*  Curve  (b)  corresponds  to  precompensation  only,  again  with 
6 * 0.9*  Curve  (c)  corre  ponds  to  tilt  removal  and  postcompensation 
(6  = 0,  near  field  propagation),  while  (d)  applies  for  postcompensation 
only.  Curve  (e)  shows  the  result  for  tilt  remova I only  (no  postcompensation), 
while  curve  (f)  corresponds  ♦:«•>  no  compensation  whatsoever.  Examination  of 
these  curves  reveals  the  gain  in  quality  factor  achieved  by  each  type  of 
image  compensation. 


Fig.  5:  Q vs.  N:  (a)  precompensation  and  post- 

compensation; (b)  precompensation  oniy; 

(c)  tilt  removal  and  postcompensation; 

(d)  postcompensat  ion  oniy;  (e)  tMt 
removal  only;  (f)  no  compensation. 
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III.  COMPARATU.,  NOISE  PERFORMANCE  OF  AMPLITUDE 

INTERFEROMETRY  AND  SPECKLE  INTERFEROMETRY 

In  this  section  we  present  an  analysis  which  compares  the  signal- 
to-noise  performance  of  two  techniques,  one  known  s ampl i tude  interfero- 
metry and  the  second  known  as  speckle  interferometry.  In  order  to  make 
this  section  as  self  contained  as  possible,  we  repeat  certain  portions 
of  material  presented  in  previous  reports  [1,2],  but  an  attempt  Is  made 
to  keep  this  repetition  to  a minimum. 

We  assume  that  we  desire  to  obtain  information  about  a distant 
object,  and  that  the  atmosphere  intervenes  between  our  observation 
instrument  (a  telescope)  and  that  object.  Both  techniques  of  interest 
are  used  to  recover  only  partial  information  about  the  object  - namely, 
we  wish  to  estimate  the  strength  of  the  squared  modulus  of  the  Fourier 
spectrum  of  the  object.  For  the  techniques  of  concern  here,  no  attempt 
is  made  to  recover  the  phase  of  the  object  spectrum. 

A.  SPECKLE  INTERFEROMETRY 

By  the  term  "speckle  interferometry",  we  mean  the  technique  intro- 
duced by  Labeyrie  [3]  which  depends  on  averaging  the  "energy  spectra" 
(squared  moduli  of  the  Fourier  transforms)  of  a sequence  of  k short- 
exposure  photographs.  By  the  term  "amplitude  interferometry",  we  mean  the 
general  class  of  techniques  which,  by  the  insertion  of  masks  or  wavefront 
folding  prisms,  operates  the  telescope  as  a Fizeau  stellar  interferometer 
or  a set  of  Fizeau  stellar  interferometers  with  different  aperture  spacings. 
Examples  are  the  interferometers  of  Currie  [^]  and  of  Breckinridge  [5]. 

Considering  first  the  case  of  speckle  interferometry,  the  predictions 
of  s ignal -to-noi se  ratio  performance  are  based  on  a specific  model  of 
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the  photodetection  process  which  takes  into  account  only  the  noise 
generated  by  the  finite  number  of  photons  utilized  in  the  recording  of 
the  image  date  and  noise  introduced  by  atmospheric  fluctuations.  The 


image  dj  (x)  detected  on  the  exposure  (exposure  time  assumed 


short  compared  with  the  fluctuation  time  of  the  atmosphere)  is  modeled  as 

,th 


an  inhomogeneous  Poisson  process.  The  form  of  the  j detected  image 
is  thus 

N . 

J . . 

(ID 


d.(x)  = l6(x-x.,) 

->  k=l 


where  each  6-function  represents  a discrete  photoevent,  and  Xj^j  is 


th 


the  location  on  the  spatial  detector  of  the  k photoevent.  N 


:th 


represents  the  total  number  of  photoevents  recorded  during  the  j 
exposure.  This  inhomogeneous  Poisson  process  is  taken  to  have  a rate 


A j (x)  which  is  proportional  to  the  classical  image  intensity 


f K 

incident  on  the  detector  during  the  j exposure  interval, 


xAt)  - iDx) 


hv  ■* 


(12) 


where 

n is  the  quantum  efficiency  of  the  detector; 

T is  the  image  measurement  time  (assumed  identical  for  all 
images) ; 

h is  Planck's  constant;  and 
V is  the  mean  optical  frequency. 

The  speckle  interferometry  technique  involves  the  following  initial 
steps : 


I.  E ch  detected  image  dj(x)  is  Fourier  transformed  to  produce 


a 2-D  spectrum  Dj (v)  , where 
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(13) 


3. 


The  squared  modulus  of  each  such  spectrum  is  taken,  producing 

I Dj ( v) I ^ ; 

These  "energy  spectra"  are  averaoed  over  the  set  of  k 
independent  pictures, 


Evaluation  of  the  expected  value  of  shows  [2]  that  it 

can  be  expressed  as 


^D(v)|^^]  = N+(N)^i.(v) 


(15) 


where 


N is  the  expected  number  of  photoevents  in  a single  detected 
image,  and 

<J),(v)  represents  the  expected  value  of  the  energy  spectrum  of 
the  classical  image  intensity  incident  on  the  detector, 
normalized  so  that  $1(0)  = 1. 

In  turn,  we  can  express  ^.(v)  as 


‘l’j(v)  = E 


O(^) 


(16) 


where  Bj  (v)  is  the  optical  transfer  function  of  the  atmosphere-telescope 
combination  on  the  exposure,  and  |0(v)|^  is  the  energy  soectrum 

of  the  true  object  radiance  distribution  o(x)  , but  normalized  such  that 

|0(0) 1^=1. 

If  we  assume  that  the  telescope  diameter  D is  much  greater  than 


_ li,  _ 


Fried's  atmospheric  coherence  diameter  r^  ')  • 

mid-frequency  region  of  the  telescope  OTF  it  is  known  [3]  that 


E 


(17) 


where  B^(v)  represents  the  telescope  OTF.  Combining  (15).  (16)  and 
(17).  we  find  that  a reasonable  unbiased  estimate  of  |0(v)|^  is 
given  by 


2 


- N 

B^(v) 


(18) 


In  assessing  the  performance  of  the  estimator  of  Eq.(l8),  we  must 

find  its  variance  for  a given  finite  k.  If  in  the  end  we  wish  to  express 

a s ignal -to-noise  ratio  associated  with  this  estimator,  we  can  equally 

well  calculate  the  s igna 1 -to-noi se  ratio  associated  with  |D(v)|^ 

Fluctuations  are  expected  in  our  estimate  due  to  the  finite  number  of 

photoevents  detected  and  due  to  fluctuations  of  the  atmosphere  from  picture 

2 

to  picture.  The  variance  a associated  with  a single  measured  value 

, ,2 

of  |D|  has  been  shown  [1,2]  to  be 

„ = N + Cn)^  + 2(2+^)  (N)^  4>.(v) 

|d|2  _ „ 

+ (N)^4>.  (2v)  + (hf)^  (v) 

2 

Removing  the  bias  from  the  expected  value  of  |D|  in  Eq.(15),  the  rms 
s ignal -to-noi se  ratio  associated  with  a single  frame  becomes 
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(20) 


1+n  B^(v) I 0(v) 


n B^(v) 1 0(v) 1 ^ 


+-^  l+4n  B_(v)!0(v)1^  + n B^(2v)j0(2v) 
N '■  ‘ 


2 


where  n = (*]^)  ^ represents  the  average  number  of  photoevents  per 

frame  per  atmospheric  coherence  area.  If  the  k detected  frames  are 
recorded  with  independent  realizations  of  the  atmospheric  transfer  function, 
the  rms  s i gna 1 - to-no i se  ratio  associated  with  k frames  becomes 


(21) 


These  results  are  consistent  with  those  obtained  eariier  by  Roddier  [6]. 

B.  AMPLITUDE  INTERFEROMETER 

For  the  case  of  amplitude  Interferometry,  several  different  measure- 
ment configurations  can  be  imagined.  Two  such  configurations  are  considered 
here,  and  the  s ignal -to-noi se  ratios  achieved  are  shown  to  be  equal. 
Nonetheless,  this  s i gna 1 -to-no i se  ratio  is  different  than  that  achieved 
by  speckle  interferometry. 

The  first  method  of  interest  is  obtained  with  only  a slight  modifi- 
cation of  the  speckle  interferometry  technique.  We  insert  a mask 
in  the  aperture  of  the  telescope;  the  mask  contains  two  circular  open- 
ings, each  of  diameter  approximately  r^  [7]  and  separated  by  a fixed 
vector  distance  rj  , as  shown  in  Fig.  6(a).  The  corresponding  diffraction- 
limited  telescope  OTF  has  the  form  shown  in  Fig.  6(b).  The  center 
frequency  Vj  of  the  two  bandpass  islands  of  the  OTF  is  given  by 
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The  bandwidth  of  these  islands  is 


B = . (23) 

A F 

and  the  peak  amplitude  is  1/2. 

A sequence  of  < pictures  is  taken,  each  with  an  exposure  time 
short  compared  to  the  fluctuation  time  of  the  atmosphere.  The  resulting 
photon-limited  images  are  Fourier  transformed,  and  the  squared  moduli 
of  these  transforms  are  averaged  over  the  set  of  pictures,  as  before. 

The  effect  of  the  atmosphere  on  any  given  picture  is  more  simple  in  this 
case  than  in  the  previous  case  of  more  general  speckle  imaging.  If  the 
scintillation  effects  introduced  by  propagation  through  the  atmosphere  are 
negligibly  small,  then,  due  to  the  fact  that  the  individual  aperture 
openings  are  of  size  r^  , the  only  effect  of  the  atmosphere  is  to 
randomly  shift  the  phase  of  all  Fourier  components  within  the  bandpass 
spectrum  by  a fixed  amount.  This  phase  shift  is  identical  for  all 
frequency  components  within  the  spectral  island,  but  it  changes  from 
frame  to  frame. 

The  signal -to-noise  ratio  associated  with  the  estimate  |0(vj)|^ 

differs  from  that  achieved  in  the  previously  considered  speckle  imaging 

case.  One  reason  for  the  difference  is  the  different  effect  that  the 

atmosphere  has  on  the  detected  images,  as  mentioned  above.  A second 

reason  arises  from  the  fact  that  the  collecting  area  of  the  speckle 

2 

imaging  system  is  tiD  /4  while  the  collecting  area  in  tlie  present  case 
2 

is  only  irr^/Z. 

If  we  remove  terms  introduced  by  the  fluctuating  amplitude  of  the 

2 

atmospherical  I y- induced  OTF,  the  variance  of  |d|  becomes 

T * 2n  + 4(n)^  + I6(n)^<t’.  + l6(n)^‘t>.  (2^) 

1D|^ 
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where  we  have  taken  account  of  the  fact  that  the  average  number  of 
photoevents  detected  per  frame  is  now  2n.  The  expected  value  of  |d| 
is  in  this  case  given  by 


2n  + 4(n)^*. 


(25) 


Removing  the  bias  from  this  expected  value,  and  noting  that 


'l'.(vj)  = |B^(vj)  l^lO(vj)  1^  = ^|0(vp|^  , (26) 

we  obtain  for  the  single-frame  rms  signal -to-noi se  ratio 


(1) 


n|0(v,)|V2 


1 + -^  + (1+n) l0(v, ) l^j 
2n  ) 


(27) 


A disadvantage  of  the  above  measurement  technique  is  that  only  one 
spacing  can  be  explored  at  one  time.  More  desirable  would  be  a 

system  that  allows  many  different  fringes  to  be  detected  at  one  time. 

Such  a system  has  been  proposed  by  Currie.  The  wavefront  from  the  full 
unobstructed  telescope  aperture  passes  through  a wavefront  folding  inter- 
ferometer, as  illustrated  in  Figure  7-  A detector  placed  at  a single 
point  at  the  output  of  the  wavefront  folding  interferometer  detects  the 


Fig.  7:  Amplitude  interferometer 
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interference  of  light  from  two  points  on  the  telescope  aperture,  those 
points  being  located  at  l80“  relative  angular  rotation  about  the  center 
of  the  telescope  mirror.  An  entire  array  of  detectors  at  the  output  of 
the  wavefront  folding  interferometer  thus  detects  simultaneously  the 
interference  between  many  different  pairs  of  points,  each  pair  being 
separated  by  a different  vector  spacing  If  a frequency  shifter  is 

placed  in  one  arm  of  the  wavefront  folding  interferometer,  each  detector 
will  see  an  oscillating  intensity,  the  amplitude  of  which  is  a direct 
measure  of  the  amplitude  of  the  object  spectrum  at  the  corresponding 
spatial  frequency. 

The  output  d^t)  of  a given  detector  is  subjected  to  an  electronic 
measurement  system  as  shown  in  Fig.  8.  The  frequency  coq  corresponds 
to  the  angular  frequency  shift  introduced  in  one  arm  of  the  interfero- 
meter. The  system  shown  detects  the  in-phase  and  quadrature  components  of 
the  oscillating  signal,  and  combines  these  to  estimate  the  squared  magnitude 
of  the  fringe  amplitude.  It  is  again  assumed  that  a single  detector  sees 

light  collected  from  two  wavefront  patches,  each  of  diameter  r . 

0 

2 2 

The  measured  quantity  Z = C + S can  be  shown  to  have  the 
following  mean  and  variance: 


Z = 2n  + (n)2|0(v)|^ 

“ 2n  + ^(rT)^  + 4(n)^|o(v)  I ^ + 4(n)^|0(v)|^ 


(28) 


Removing  the  b'as  from  Z , we  obtain  an  rms  signal -to-noise  ratio  (based 
on  a single  integration  interval  that  is  short  compared  to  the  fluctuation 
time  of  the  atmosphere) 
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Fig.  8:  Fringe  contrast  detector 


Z - 2n  ^ ni0(v)|^/2 

^ ""Z  j,  + i_+  (,.^-)|3(^)|2j^ 

( 2n  ) 

which  is  identical  with  the  result  obtained  previously. 
C.  COMPARISON  OF  SPECKLE  AND  AMPLITUDE  INTERFEROMETRY 


(29) 


We  now  compare  the  single-frame  rms  s ignal -to-noi se  ratios  appropri- 
ate for  speckle  and  amplitude  interferometry.  In  making  the  comparison, 
we  assume  that  the  optical  bandwidths,  integration  times  and  optical 
efficiencies  of  the  two  systems  are  identical.  In  addition,  we  assume  [8] 
that  the  detectors  for  the  two  systems  have  identical  quantum  efficiencies. 

Figure  9 shows  plots  of  the  single-frame  rms  s i gna 1 -to-no i se  ratio 
for  a point-source  object  (|0(v)l^=  I).  For  the  case  of  the  speckle 
interferometer  we  assume  that  v/v^  = 0.8  (vq  - diffraction-limited 
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Fig.  9:  Signal-to-noise  ratios,  speckle 

vs.  amplitude  interferometer, 

|6|2  = 1 


cutoff  frequency)  and  (1.6  meter  telescope  d’ameter  with 

r^  = 10  cm).  The  chief  conclusions  to  be  drawn  from  these  curves  are: 

1.  For  n <<  1 , the  rms  signal-to-noise  ratios  of  both  techniques 

7/2 

increase  in  proportion  to  (n) 

sO) 

ampl i tude  — 

2/2 

,p«kle- 

N rg  T 

For  the  parameters  chosen  here  speckle  interferometry  is  superior 
to  ampiitude  interferometry  in  this  range  by  a factor  of  4.7; 


'1 


2.  For  n in  the  range  O.i  < n < 1.0  , the  curves  cross,  and 
amplitude  interferometry  becomes  superior  to  speckle 

i nterferometry ; 

3.  For  n >>  1 , the  •j;j-  for  speckle  interferometry  asymptotically 
approaches  unity,  never  exceeding  that  value,  while  for  amplitude 
interferometry  it  continues  to  increase  as  n/TT  h i ndef i ni tely . 

In  Figure  10  we  show  the  corresponding  curves  for  “ 0.8, 

D/r^  = 16  , but  a much  smaller  object  spectrum  |o|  = 10  . For  very 

small  values  of  n , speckle  interferometry  remains  superior  to  amplitude 

interferometry,  but  by  a smaller  margin  than  previously.  In  the  vicinity 

— -2  — 
of  n = 10  for  speckle  interferometry  and  n = 1 for  amplitude  inter- 


Fig.  10;  S i gnal -to-noi se  ratio,  speckle 
vs. 

Ol*^  - 10- 


amplitude  interferometry 

I 2 _ , n-H 
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ferometry,  both  curves  change  from  an  (n)'^  increase  to  an  n 
increase,  with  amplitude  interferometry  superior  to  speckle  interfero- 
metry. Eventually,  for  values  of  n so  large  as  to  be  off  this  figure, 
the  speckle  interferometry  curve  will  saturate  at  value  unity,  while  the 
amplitude  interferometry  curve  will  begin  to  increase  in  proportion  to 


IV.  SPACE-VARIANT  AND  NONLINEAR  RESTORATION  TECHNIQUES 

Our  past  work  on  the  restoration  of  photon-limited,  blurred  Images 
has  considered  only  the  case  of  linear,  space- i nvar iant  restoration 
filters.  Based  on  the  intuitive  feeling  that  improved  restorations  can 
be  obtained  with  linear  space-variant  filters  and  nonlinear  filters,  we 
have  begun  to  explore  such  possibilities.  The  question  which  is  of 
major  concern  here  is  how  much  improvement  in  restoration  filter  perfo 
ance  can  be  achieved  by  going  to  the  decidedly  more  complex  space- 
variant  and  nonlinear  restoration  schemes.  Our  work  in  this  area  has 
progressed  at  least  to  the  point  where  the  formulation  of  the  problem  has 
become  clearer. 

In  anticipation  that  future  technical  reports  will  detail  this  work 
more  thoroughly,  we  mention  here  the  current  state  of  our  knowledge, 
without  any  specific  "proofs"  of  the  results  stated. 

The  most  general  formulation  of  the  linear  least-squares  filtering 
problem  can  be  appreciated  by  reference  to  Figure  11.  Foi  simplicity, 
we  treat  space  as  a one-dimensional  variable  x , the  results  being 
easily  generalized  to  two  dimensions.  The  object  radiance  distribution 
o(x)  is  assumed  to  b i a random  process,  but  for  generality  it  is  allowed 
to  have  space-variant  or  non-stat ionary  statistics.  The  object  o(x) 
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o(x) 
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Fig.  II:  Space-variant  least-squares 

f i 1 ter ing 

is  degraded  by  a linear  blur  with  space-variant  point  spread  function 
b(x,x"‘).  This  blurred  image  is  then  converted  to  a detected  image  by 
passage  through  the  "Poisson  generator",  which  generates  an  inhomogeneous 
Poisson  impulse  process,  with  rate  A(x)  proportional  to  the  classical 
image  intensity, 

A(x)  - ^\{x)  (30) 

h V 

The  detected  image  then  passes  through  a linear,  possibly  space-variant, 
restoration  filter  with  impulse  response  h(x,x").  The  impulse  response 
h(x,x')  is  chosen  to  minimize  the  mean  square  error  between  the  restored 
image  r(x)  and  an  "ideally  filtereu  ' version  of  the  object,  o(x) . The 
ideal  filter  may  in  general  be  space-variant,  and  its  impulse  response  is 
represented  by  s(x,x"). 
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Without  presenting  a proof,  we  simply  state  that  the  impulse 
response  of  the  optimum  restoration  filter  can  be  shown  to  he  a function 


h(x,x")  satisfying  the  integral  equation 


h(x,x'')R^(x",x"')dx'" 


(31) 


s(x,x")R^^(x',x")dx" 


where 


R^(x",x"')  = E[d(x")d(x")  ] 


(32) 


= E[d(x')o(x"")]  . 


In  turn,  the  correlation  functions  of  interest  can  be  expressed  as 


Rj(x'.x'-') 


2^2 


i (x)6(x'-x")  + R.  (x",x'0 

h'  h V ' 


(33) 


Rdo(x%x^') 


hv 


where 


R.(x",x'")  = E[i (x") i (x"")J 


(34) 


R.q(x'.x'^) 


E[i  (x')o(x")] 


Finally  we  note  that 


i(x) 


b(x,x'')o(x'")dx"  . 


(35) 


A case  of  major  interest  is  obtained  if  the  blur  spread  function 
b and  the  ideal  spread  function  s are  space- invar  iant , but  the  object 
itself  is  non-stat ionary . Since  the  objects  of  prime  interest  here  are 
of  finite  size,  they  must  of  necessity  have  spatially  non-stat ionary 
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statistics.  The  correlation  functions  of  Interest  become 


R.  (x"  ,x") 


b(x'-5)b{x""-n)RQ(C,n)d5dn 


R.  (x',x"‘) 
10 


b(x'-?)RQ(5,x-)d& 


(36) 


where  R (C,x"')  is  the  autocorrelation  function  of  the  non-stat ionary 
o 

object . 

In  some  cases  it  is  possible  to  argue  that  the  object  is  non-stat ionary 
primarily  by  virtue  of  a non-stat ionary  second  moment.  In  such  cases, 
the  autocorrelation  function  of  the  object  takes  the  approximate  form 


(37) 


The  idealized  case  of  an  object  with  a "white"  spectrum  yields  the  auto- 
correlation function 


R^(x',x^O  s 6(x--x-)  . (38) 

An  additional  case  of  interest  is  that  of  a stationaty  object 
multiplied  by  a "window"  function  w(x).  If  r^(x'-x^^)  represents  the 
autocorrelation  function  of  the  original  stationary  object,  the  auto- 
correlation of  the  windowed  object  is 

R^(x",x")  = w(x")w(x"')r^(x-'-x")  . (39) 

It  is  perhaps  worth  special  mention  here  that,  if  the  blur  filter 
and  the  ideal  filter  are  space- invar i ant  and  the  object  i s stat i st ical 1 y 
stationary,  then  we  can  show  from  the  results  presented  above  that  the 
filter  which  achieves  minimum  mean-square  error  is  always  the  space- 
invariant  filter  derived  in  our  previous  reports  [1,2].  Such  a filter 
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remains  optimum  in  spite  of  the  s ignal -dependent  character  of  Poisson 
no  I e. 

In  addition  to  the  linear  space-variant  filters  discussed  above, 
we  are  also  beginning  preliminary  studies  of  the  nonlinear  filter  shown 
in  Fig.  12(a).  In  this  case  the  filter  is  taken  to  be  linear  and  space 
invariant,  with  the  transfer  function  (see  Eq.(4)) 


Fig.  12:  Realization  of  non-linear  filter 

(a)  direct;  (b)  with  feedback 

- 28  - 


H{^) 


kN  S B * 4>q 
i+n|b|^$q 


appropriate  for  minimizing  mean-squared  error.  However,  the  parameter 
N i s no  longer  taken  to  be  equal  to  the  average  number  of  photoevents 
in  the  entire  image.  Rather,  the  detected  image  d(x)  is  passed  through 
a linear,  invariant  smoothing  filter  with  impulse  response  g(x)  , the 
purpose  of  which  is  to  produce  a local  average  value  of  N.  Thus  N 
is  a function  of  x , and  the  filter  parameter  is  chosen  appropriately 
for  the  local  photon  environment.  In  regions  of  the  picture  where  there 
have  been  very  few  photoevents,  very  little  enhancement  will  be  attempted, 
but  in  regions  where  the  number  of  photoevents  is  large,  greater  enhance- 
ment is  attempted.  Such  a filter  has  a characteristic  that  is  highly 
dependent  on  the  detected  signal  d(x)  ; it  is  this  signal -dependent 
character  that  causes  the  filter  to  be  nonlinear.  There  is  no  doubt  that 
such  a filter  can  perform  better  than  any  linear  filter  in  restoring  photon- 
limited  images.  The  key  questions  are  how  much  better  it  can  perform, 
and  whether  the  increased  performance  is  worth  the  extra  computational 
complexity  involved.  As  a potentially  helpful  fact,  we  note  that  the 
required  transfer  function  can  be  realized  by  the  feedback  structure  of 
Fig.  12(b),  where  the  parameter  N is  controlled  simply  by  varying  the 
gain  of  the  feed-forward  filter. 

Work  on  the  topics  outlined  in  this  section  is  still  in  progress  and 
will  be  reported  on  in  greater  detail  in  our  next  technical  report. 
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MmUC  SYSTEM 


■A8E  UNITS: 

Quantity  UnM  SI  Symbol  . Ponwuto 


length 

metre 

m 

... 

maw 

kilogram 

time 

aecond 

• 

electric  currant 

ampere 

A 

thetmodynamlc  temperature 

kelvin 

K 

amount  of  lubetanca 

mole 

owl 

lumlnoue  Intanalty 

candela 

cd 

SUrPLEMENTARY  UNITS: 

plana  angle 

radian 

rad 

lolld  angle 

ataradian 

ir 

... 

DERIVED  UNITS: 

Accalaration 

metre  per  aecond  aquarad 

aVa 

activity  (of  a radioactive  aouice) 

diaintegratlon  par  aecond 

(dlaintagrallon|/a 

angular  accalaration 

radian  par  aecond  aquared 

nd/a 

angular  velocity 

radian  per  aecond 

lad/a 

area 

aquare  metre 

m 

denaity 

kilogram  per  cubic  metre 

kg/m 

electric  capacitance 

farad 

F 

AuA' 

electrical  conductance 

aiemena 

S 

AA/ 

electric  flald  etrength 

volt  per  metre 

V/m 

electric  Inductance 

henry 

H 

V.a/A 

electric  potential  different 

volt 

V 

W/A 

electric  reslitance 

ohm 

V/A 

electromotive  force 

volt 

V 

W/A 

energy 

ioule 

1 

N.m 

entropy 

loule  per  kelvin 

... 

VK 

force 

newton 

N 

kg>m/a 

frequency 

hertz 

Hz 

(cyciaVa 

illuminance 

lux 

lx 

Im/m 

luminance 

candela  par  aquare  metre 

... 

od/m 

lumlnoue  flux 

lumen 

Im 

cd.fr 

magnetic  field  atrength 

ampere  par  metre 

... 

A/m 

magnetic  flux 

waber 

Wb 

V-a 

magnetic  flux  denaity 

taala 

T 

Wbfrn 

magiMtomotive  force 

ampere 

A 

... 

power 

watt 

W 

V* 

preaaura 

paacal 

Pa 

Him 

quantity  of  electricity 

coulomb 

C 

A*a 

quantity  of  heat 

ioule 

) 

N.m 

radiant  intenaity 

watt  per  ateradian 

W/ar 

apecific  heat 

ioule  per  kilogram-kelvin 

i/kg.K 

atreaa 

paacal 

Pa 

N/m 

thermal  conductivity 

wat‘.  per  metra-kalvin 

W/tn-K 

velocity 

metre  par  aecond 

m/a 

viacoaity,  dynamic 

paacal-aacond 

Pa-a 

vlacoalty,  kinematic 

aquare  metre  per  aecond 

m/a 

voltage 

volt 

V 

W/A 

volume 

cubic  metre 

m 

wavenumber 

reciprocal  metre 

(wave)/m 

work 

ioule 

i 

N-m 

SI  PREFIXES: 


Multiplication  Faclort 


Prafix  SI  Syinbul 


1 000  000  000  000  > 10“ 

1 000  000  000  • 10* 
1000  000-  10* 

1 000-  10> 
100  - 10< 
10  - 10> 
0.1  - 10-' 
0.01  - 10-' 
0.001  - 10-» 
0.000  001  - 10-* 
0.000  000  001  - 10-* 
0.000  000  000  001  - 10-“ 
0.000  000  000  000  001  - 10-“ 
0.000  000  000  000  000  001  - 10-“ 
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MISSION 

of 

Rom  Air  Development  Center 


RADC  plans  and  conducts  research,  exploratory  and  advanced 
development  programs  in  command,  control,  and  coummications 
(C^)  activities , and  in  the  areas  of  information  sciences 
and  intelligence . The  principal  technical  mission  areas 
are  communications , electromagnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microeave 
physics  and  electronic  reliability , maintainability  and 
cong>atibility . 


